clc; close('all'); clear;
load('Datos11164');
T = 0.1;
tn = 50;        % Segundos
N = tn / T;     % Numero de muestras para tn Segundos
time = a_ent_esc(:,1);
in = a_ent_esc(:,2);
out = b_ent_esc(:,2);
txtB = {'Entrada','Respuesta en C','Recta',...
    'Estimacion del Modelo G1(s)'};
[Tu, Ta, K] = identificar(time, in, out, N, txtB);
T1 = Ta / 4;
T2 = 2*T1;
s = tf('s');
G1 = K /((1+T1*s)*(1+T2*s));
